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The  cost  and  availability  of  platinum,  together  with  degradation  phenomena  that  are  directly  or  indirectly 
caused  by  the  use  of  a  precious-metal  catalyst,  are  major  obstacles  to  the  commercialization  of  incumbent 
polymer  electrolyte  membrane  fuel  cell  (PEMFC)  technologies.  The  Flowcath®  fuel  cell  technology  is  a 
highly  promising  alternative  to  the  conventional  PEMFC  that  avoids  many  of  these  issues.  While  the 
anode  resembles  that  of  a  conventional  cell,  the  novel  cathode  employs  a  liquid-based  polyoxometalate 
redox  mediator  and  is  entirely  free  of  any  precious  metals.  A  transition-metal  based  redox  centre  within 
the  polyoxometalate  undergoes  reversible  reduction  in  the  cathode  and  is  regenerated  externally  using 
air.  In  this  paper,  the  performance  of  the  cell  for  different  catholyte  flow  rates  and  polyoxometalate 
concentrations  is  described  and  a  first  model  is  developed  and  validated.  The  voltage  losses  are  quantified 
and  compared  to  the  values  for  a  conventional  cell  employing  a  standard  membrane  electrode  assembly. 
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1.  Introduction 

There  are  several  impediments  to  realising  the  commercial  ben¬ 
efits  of  proton  exchange  membrane  (PEM)  fuel  cells,  notably  cost 
and  longevity.  The  use  of  platinum  in  conventional  PEM  fuel  cells, 
particularly  in  the  performance-limiting  cathode,  not  only  raises 
the  cost  of  the  cells  but  leads  directly  to  degradation  phenomena 
such  as  dissolution  and  sintering  of  the  catalyst  particles.  A  PEM 
fuel  cell  that  replaces  the  conventional  cathode  with  a  flow-through 
platinum-free  carbon  electrode  was  developed  by  ACAL  Energy  Ltd. 
(Fuel  Cells,  Patent:  WO/2007/1 10663).  This  novel  structure  utilises 
a  polyoxometalate  (POM)  containing  a  transition  metal  ion  redox 
centre,  such  as  Viv/Vv,  in  place  of  hexavalent  tungsten  or  molybde¬ 
num.  The  basic  structure  of  the  POM  anions  used  in  the  Flowcath® 
system  is  of  type  [XMi2_mVm04o]n_,  where  X  is  commonly  P,  Si  or 
B,  M  is  typically  molybdenum  or  tungsten,  and  n  is  the  charge  of 
the  anion.  The  fundamental  reaction  in  the  cathode  is: 

(XM12_m<O40)n-  +xe-  +xH+  — ►  (HxXM12_mVV_XvO40)n- 

(1) 
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where  Viv  and  Vv  are  the  reduced  and  oxidised  (vanadium) 
species,  respectively,  and  x  is  between  1  and  m  depending  on  the 
extent  of  the  reaction.  In  a  separate  regenerator,  the  POM  is  oxi¬ 
dised  in  the  solution  phase  before  it  is  recycled  to  the  electrode: 
(HxXMu_mV"(04o)n-  +(x/4)02  (XMU-mVmO40)n-  +(x/2)H20. 

In  this  paper,  the  performance  of  the  system  is  presented  for  the 
first  time.  To  aid  understanding  of  the  performance,  a  first  model 
of  the  cell  (focused  on  the  cathode  processes)  is  developed  and  the 
cell  voltage  losses  for  different  flow  rates  and  POM  concentrations 
are  studied. 

2.  Experimental  and  modelling  details 

2.1.  Experimental 

A  single  25  cm2  cell  was  tested  with  the  standard  ACAL  catholyte 
(POM  and  de-ionised  water)  at  varying  flow  rates  and  catholyte 
concentrations.  The  cell  build  incorporated  a  porous  carbon  cath¬ 
ode  (RVC  supplied  by  ERG)  material  for  the  reduction  of  the 
catholyte,  a  commercially  available  50  pirn  Nation  based  membrane 
(supplied  by  Ion  Power)  and  a  gas  diffusion  layer  (GDL)  with  a 
flow  field  and  microporous  layer  at  the  anode.  These  were  encap¬ 
sulated  by  two  carbon  blocks  (Le  Carbone)  and  current  collectors 
(gold  coated  copper  plates),  and  held  together  with  two  steel  blocks 
bolted  to  a  set  torque.  The  catholyte  solution  was  heated  to  80  °C 
and  pumped  through  the  cell  (KNF  diaphragm  pump)  where  the 
temperature  was  maintained  with  a  heater  rod.  H2  (99.999%  purity, 
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BOC)  was  supplied  to  the  anode  and  a  back  pressure  of  1 .6  bar  main¬ 
tained  throughout.  Electrochemical  measurements  were  taken  by 
connecting  the  cell  to  a  potentiostat  (Biologic)  and  carrying  out 
polarisation  experiments  at  500  mAs-1  to  evaluate  the  cell  volt¬ 
age.  The  polarisation  curves  were  taken  (with  the  catholyte  in  the 
same  oxidation  state)  at  50,  100  and  150  ml  min-1  for  POM  con¬ 
centrations  of  0.45, 0.3  and  0.2  M.  The  catholyte  concentration  was 
reduced  by  the  addition  of  water.  Excess  water  was  removed  by 
evaporation  in  order  to  maintain  a  constant  concentration  during 
testing. 


2.2.  Mathematical  model 


Since  the  conventional  anode  is  well-characterised,  attention  is 
focused  on  the  complex  cathode  processes.  The  anode  is  treated 
implicitly  through  boundary  conditions  and,  along  with  contact 
resistances,  as  an  additional  resistance  Ra  that  gives  rise  to  a  cell 
voltage  loss  IRa ,  where  /  is  the  applied  cell  current  (load).  This  is  a 
reasonable  first  approximation  since  the  activation  and  mass  trans¬ 
port  losses  in  a  conventional  anode  are  small  and  often  assumed  to 
be  zero.  The  domain  to  be  modelled  is  enclosed  in  the  dotted  rect¬ 
angle  in  Fig.  1.  Based  on  the  POM  used  in  the  experiments,  reaction 
(1)  is  represented  in  the  following  manner: 

Vv  +  e_  — >  Viv  (2) 


where  Viv  and  Vv  refer  to  the  redox  centres  of  the  POM.  The 
species  considered  are  dissolved  H20,  H+,  Viv,  Vv  and  a  counte¬ 
rion.  Certain  simplifying,  (but  reasonable)  first-approximations  are 
adopted:  (i)  two-dimensional  through-plane  flow  (see  Fig.  1);  (ii) 
a  dilute-solution  approximation  [1];  (iii)  laminar,  incompressible 
flow;  (iv)  cathodic  and  anodic  charge  transfer  coefficients  of  0.5; 
(v)  the  POM  solution  is  returned  to  the  electrode  from  the  regener¬ 
ator  at  a  constant  POM  concentration;  (vi)  reduction  of  vanadium 
species  within  the  POM  molecule  is  a  one  electron  process  (for  a 
similar  system  Odyakov  et  al.  [4]  demonstrated  that  this  assump¬ 
tion  is  valid). 

Using  the  Nernst-Planck  equation  [1],  the  mass  balances  for 
H20,  Vv,  Viv  and  H+  in  the  porous  carbon  cathode  take  the 
form: 


v- 


jh- C3/2D/Vq 


ZiCi€3/2Di 

RT 


FV0e  -  Cj 


Vjjc 

F 


(3) 


where  6  is  the  porosity;  0e  is  the  ionic  potential;  R,  T and  F  are  the 
molar  gas  constant,  temperature  and  Faraday’s  constant,  respec¬ 
tively;  Cj,  D/,  Zj  and  v*  are  the  concentration,  diffusion  coefficient 
(subject  to  a  Bruggeman  correction),  charge  and  stoichiometric 
coefficient  (in  reaction  (2))  of  species  i;  and  jc  is  the  cathode  transfer 
current  density.  The  counterion  concentration  is  given  by  the  elec¬ 
troneutrality  condition:  J2izici  =  0.  The  electrolyte  velocity  v  (in  the 
last  term  on  the  left-hand  side)  is  given  by  Darcy’s  law  in  combina¬ 
tion  with  the  Kozeny-Carman  law  for  the  electrode  permeability 

[2] ,  where  /x  is  the  dynamic  viscosity  of  the  liquid,  dp  is  the  mean 
pore  diameter  and  I<  is  the  Kozeny-Carman  constant.  Ionic  charge 
conservation  in  the  cathode  electrolyte  can  be  expressed  as  follows 

[3] : 


jc  =  -V  •  \^V<PeY/zfD,£y2Ci  +  F^/ziDf  Vqj  (4) 

and  the  electronic  potential  </>s  in  the  porous  electrode  and  the 
current  collector  is  given  by  Ohm’s  law: 


-orseffV2</>s  =  -jc,  -crcV2<t>s  =  0 


(5) 


where  af  (<rc)  is  the  effective  conductivity  of  the  cathode  (collec- 
tor). 

The  driving  forces  for  the  movement  of  dissolved  protonated 
water  in  the  membrane  are  concentration,  potential  and  pressure 
gradients.  The  concentration  Ch2o  satisfies  the  mass  balance[5]: 

— V  •  [D^J Vcm2o)  +  V  ■  (vCh2o)  =  0,  i )=-lStFcH+V<pe-l^Vp  (6) 

H' 

in  which  De$  is  the  effective  water  diffusion  coefficient,  v  is  the  liq¬ 
uid  velocity  (given  by  Schloegl’s  equation),  is  the  electrokinetic 
permeability  and  kp  is  the  hydraulic  permeability.  Incompress¬ 
ibility,  i.e.,  V  •  v  =  0,  is  assumed.  Charge  conservation  is  again 
given  by  Ohm’s  law  in  the  form  [5]:  -om  V20e  =  O,  where  om  = 
F2D™+ch+  /{RT)  is  the  membrane  conductivity. 

The  rate  of  reaction  (2)  is  assumed  to  follow  the  Butler-Volmer 
law  (with  equal  transfer  coefficients): 

i-ffKc>)>'Vv,)',2{eW(|L)  -«xp(-||)}  (7) 

where  k  is  the  reaction  rate  constant  and  r\  =  0S  -  </>e  -  F°  is  the 
overpotential  in  the  cathode.  The  reversible  open  circuit  elec¬ 
trode  potential  is  estimated  using  the  Nernst  equation:  E°  =  E0'  + 
( RT / F )  In  {csvV/csyIV ),  in  which  the  effect  of  pH  variations  is  neglected. 
The  quantities  c?  are  the  concentrations  of  the  vanadium  species 
at  the  liquid-solid  interfaces  in  the  porous  cathode.  They  can  be 
related  to  the  bulk  values,  q,  by  approximately  balancing  the  rate 
of  reaction  with  the  rate  of  diffusion  of  reactant  to  the  electrode 
surface  at  steady  state  [3]: 

cs  =  Cj  +  eke-F(^-^-£0)/(2RT)(cv[V/yvv  +  cvv/yvlv) 

C''  “  1  +  ek  (e-f(<fe-&-E°)/(2RT)/KvV  +  eF(^-0e-£°)/(2RT)/Kv|v)  j 

where  i  is  Vv  or  VIv.  At  the  reaction  surface,  a  stagnant  boundary 
layer  is  incorporated.  Following  [6],  the  thickness  of  the  boundary 
layer,  8 ,  is  inversely  proportional  to  the  square  root  of  the  local  elec¬ 
trolyte  speed,  \v\,  i.e.,  8  =  k$/y/$\,  which  gives  y2-  =  D2/(/<5/ ^/\^\), 
where  k8  is  a  constant  of  proportionality. 

Fig.  1  shows  the  boundary  conditions.  At  the  cathode  inlet, 
the  liquid  electrolyte  enters  with  a  prescribed  bulk  velocity.  At 
the  outlet,  a  fully-developed  flow  condition  is  used.  The  cur¬ 
rent,  I  enters  uniformly  through  current  collectors  on  the  cathode 
side.  The  complementary  boundary  condition  at  x  =  Xo  is  derived 
from  a  one-dimensional  charge  balance  across  the  anode  catalyst 
layer.  Conservation  of  charge  demands  that  ae  dx4>e  +  crs  dxcj)s  =  Ila, 
where  as  and  <re  are  the  electronic  and  ionic  conductivities  in  the 
anode  catalyst  layer,  respectively.  Applying  this  equation  to  the 
membrane/catalyst  layer  interface  (zero  electron  flux)  gives  the 
condition  at  d£2\  in  Fig.  1.  The  ionic  conductivity  cre  is  approxi¬ 
mated  as  <7e  =  €^2F2Dh+ch+  where  €e  is  the  volume  fraction  of 
ionomer  in  the  anode. 

3.  Results  and  discussion 

The  system  of  equations  and  boundary  conditions  was  solved 
using  COMSOL  Multiphysics.  Due  to  the  lack  of  available  data  for 
this  new  system,  it  was  necessary  to  estimate  certain  parameters 
in  the  model:  the  reaction  constant,  /<,  the  proportionality  constant 
for  the  diffusion  boundary  layer  thickness,  k8,  and  the  resistance  Ra. 
The  parameter  estimation  (implemented  in  Matlab)  was  performed 
rigorously  through  a  nonlinear  least-squares  fit  of  six  simulated 
polarisation  curves  for  different  POM  concentrations  and  flow  rates 
to  experimental  data  obtained  under  the  same  conditions.  Six  data 
points  from  each  of  the  six  curves  were  used  to  form  the  least- 
squares  sum.  The  resulting  values  are  given  in  Table  1,  along  with 
the  other  parameter  values.  The  inlet  concentrations  of  Viv,  Vv 


R.  Singh  et  al.  /  Journal  of  Power  Sources  201  (2012)  159-163 


161 


Boundary  conditions 


30,:  <7(,V (j)e  ■  n  =  I/a;  Vp  n  =  0  3Q7: 

30,:  Reactant/electron  insulation  Vc.  •  n  =  0;  p  =  p 


30  3:  Reactant/proton  insulation  30  s: 


Vpn-0 

30  4:  -<7cVft  •«  =  //« 

30  5  and  30  6:  Reactant  insulation 
n  =  0 


c,=cf;  v  =  (0,vj 
30  g  and  30  l0: 

-aV0s-h=  I /a 


Fig.  1.  A  schematic  of  the  Flowcath®  cell  and  of  the  components  modelled  in  this  report  (enclosed  in  the  dotted  line). 


(strictly  speaking  of  the  redox  centres),  H+  and  H20  were  calcu¬ 
lated  based  on  the  POM  concentration,  the  pH  of  the  solution, 
and  the  degree  of  reduction  (DOR),  defined  as  cviv /( cviv  +  cvv).  For 
example,  for  a  /3  M  POM  concentration  and  a  DOR  of  d,  the  Viv  con¬ 
centration  is  fidm  M,  where  m  is  the  number  of  V  redox  centres.  For 
reasons  of  proprietary,  the  actual  Viv  and  Vv  concentrations  are  not 
shown. 

Fig.  2  compares  the  experimental  and  simulated  polarisation 
curves  at  different  electrolyte  inlet  flow  rates  with  a  fixed  POM  con¬ 
centration  and  at  different  POM  concentrations  with  a  fixed  inlet 
flow  rate.  The  model  captures  the  variations  in  the  cell  voltage  with 
changes  in  the  POM  concentration  and  flow  rate  well,  both  qual¬ 
itatively  and  quantitatively.  The  simulated  polarisation  curves  in 
Fig.  2(b)  overpredict  the  cell  voltage  at  low  current  densities,  par¬ 
ticularly  as  the  POM  concentration  is  increased.  This  is  likely  to  be 
due  to  the  assumption  that  the  reaction  rate  and  Nernst  potential 
are  independent  of  the  pH  (the  precise  stoichiometry  of  the  cath¬ 
ode  reaction  is  still  under  study).  Analysis  of  a  similar  POM  system 
[4]  has  shown  that  the  pH  of  the  POM  solution  decreases  with  an 
increase  in  the  POM  concentration  for  a  fixed  DOR.  In  other  words, 
the  stoichiometry  of  the  cathode  reaction  with  respect  to  protons 
depends  in  a  complex  manner  on  both  the  POM  concentration  and 
the  DOR;  the  present  model  does  not  incorporate  this  dependence. 


The  open  circuit  voltage  for  the  cell  (Fig.  2)  is  approximately 
0.83  V,  which  is  around  0.1  V  lower  than  that  of  a  PEM  fuel  using 
a  standard  membrane  electrode  assembly  (MEA).  The  theoreti¬ 
cal  open-circuit  potential  for  a  conventional  cell,  around  1.23  V  at 
298  K,  is  not  attained  due  to  a  mixed  potential  at  the  cathode  and, 
to  a  lesser  extent,  H2  crossover  [12].  For  a  0.45  M  POM  solution  and 
an  inlet  flow  rate  of  150  ml  min-1,  the  total  voltage  losses  at  0.7 
and  1  A  cm-2  are  0.279  V  and  0.359  V,  respectively.  The  equivalent 
numbers  for  a  standard  cell  (also  with  a  50  p,m  PFSA  membrane 
and  a  25  cm2  active  area)  can  be  considerably  higher,  particularly 
at  moderate  to  high  current  densities.  For  example,  a  cell  operat¬ 
ing  at  70  °C  on  fully  humidified  air  (stoichiometric  coefficient  of 
5)  and  H2  at  0%  RH  (as  in  the  present  system)  with  a  platinum 
loading  of  0.57  mg  cm-2  [7]  exhibited  voltage  losses  of  0.400  V  and 
0.583  V  at  0.7  and  1 A  cm-2,  respectively.  The  best  performance  was 
attained  by  humidifying  both  channels  at  80%  RH,  leading  to  voltage 
losses  of  0.321  V  at  0.7  A  cm-2  and  0.457  V  at  1  A  cm-2  -  still  signif¬ 
icantly  higher.  On  the  other  hand,  a  state-of-the-art  cell  employed 
by  Gasteiger  et  al.  [13]  with  a  25  pun  PFSA  membrane  and  a  50  cm2 
active  area  (platinum  loading  of  0.4  mg  cm-2)  exhibited  voltage 
drops  of  0.24  V  and  0.3  V  at  the  same  current  densities.  The  power 
densities  at  0.7,  1  and  1.2  A  cm-2  for  the  conventional  cell  in  [7] 
were  0.445,  0.488  and  0.171  W cm-2,  respectively.  In  the  present 
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Table  1 

Default  parameter  values  for  the  simulations  (inlet  concentrations  are  based  on  a 
0.3  M  POM  solution). 


Symbol 

Quantity 

Size 

6e 

Ionomer  volume  fraction  in  anode 

0.35 

€ 

Porosity  of  the  cathode 

0.85 

DyW 

Viv  diffusion  coefficient  in 
electrolyte 

S^xIO-^itPs-1  [8] 

DyV 

Vv  diffusion  coefficient  in 
electrolyte 

5.65  x  10-10  m2  s— 1  [8] 

Dh2o 

Water  diffusion  coefficient  in 
electrolyte 

2.3  x  10-9  m2  s-1  [9] 

Dw 

Water  diffusion  coefficient  in  the 
membrane 

5.75  x  10-10  m2  s-1  [10] 

Dn+ 

Proton  diffusion  coefficient  in  the 
membrane 

1.4  x  10-9  m2  s-1  [5] 

Dh+ 

Proton  diffusion  coefficient  in  the 
electrolyte 

9.7  x  10-9  m2  s-1  [11] 

crc 

Current  collector  conductivity 

8.33  x  104Sm-1 

creff 

us 

Porous  cathode  conductivity 

600  Sm-1 

I< 

Kozeny-Carman  constant:  porous 
electrode 

5.55 

K(t> 

Electrokinetic  permeability: 
membrane 

1.13  x  10-19  m2  [5] 

Kp 

Hydraulic  permeability:  membrane 

1.58  x  10-18  m2  [5] 

k 

Rate  constant  for  reaction  (2)a 

4.7  x  10“3  S”1 

ks 

Diffusion  boundary  layer  thickness 
constant3 

3.67  xlO-7  m3/2  s-1 

Ra 

Anode  and  contact  resistance3 

0.0069  Q 

Eo' 

Formal  potential:  Viv/Vv 

1.004  V  vs.  NHE 

cf 

Fixed  charge  site  (sulfonate) 
concentration 

1800  mol  m“3 

c® 

h2o 

Water  concentration  at  cathode 
inlet 

4.2  x  103  mol  m-3 

c° 

H+ 

Initial  H+  concentration 

500  mol  m-3 

T 

Cell  temperature 

80°C 

P  out 

Cathode  outlet  pressure 

lOOkPa 

a  Estimate  based  on  a  nonlinear  least-squares  analysis. 


system,  the  equivalent  values  for  a  0.45  M  POM  concentration 
and  150  ml  min-1  flow  rate  were  0.378,  0.461  and  0.493  W  cm-2, 
demonstrating  superior  performance  at  high  current  densities  (typ¬ 
ical  of  automobile  applications)  and  comparable  performance  at 
moderate  current  densities.  The  state-of-the-art  system  in  [13] 
achieved  power  density  values  of  0.50,  0.66  and  0.76  W  cm-2  and 
did  not  suffer  severe  mass-transport  losses  up  to  1.2  A  cm-2  (see 
below).  Improvements  are  clearly  required  in  order  to  match  this 
level  of  performance.  It  is  important  to  bear  in  mind,  however,  that 
the  performance  of  the  test  cell  used  in  this  work  is  far  from  optimal 
and  that  significant  improvements  are  achievable,  particularly  with 
respect  to  the  ohmic  losses  incurred  in  the  membrane,  electrodes 
and  current  collectors,  the  flow  distribution  through  the  electrodes, 
contact  resistances  and  the  POM  solution  composition. 

Fig.  2(a)  shows  that  the  performance  improves  dramatically 
with  an  increase  in  the  flow  rate  from  50  to  100  ml  min-1  but 
only  mildly  with  an  increase  from  100  to  150  ml  min-1.  The  flow 
rate  influences  both  the  characteristic  residence  time  for  reaction, 
r R  =  h /vy  (where  h  is  the  electrode  height  and  vy  is  the  y  component 
of  the  velocity)  and  the  characteristic  diffusion  time  through  the 
boundary  layer,  rD  =  82/nDi  =  k2/nvyDi ,  where  n  is  a  constant  that 
depends  on  the  pore  geometry.  The  timescale  rR  increases  rapidly 
as  the  flow  rate  is  reduced,  which  leads  to  greater  depletion  of  the 
Vv  reactant  at  a  fixed  current  density  and,  therefore,  higher  average 
cathode  overpotentials.  This  effect  is  compounded  by  the  simul¬ 
taneous  increase  in  rDt  which  leads  to  a  greater  deviation  of  the 
surface  concentration  from  the  bulk  value.  The  derivatives  of  both 
time  scales  with  respect  to  vy  are  proportional  to  1  / v y,  suggesting 
that  improvements  in  performance  decay  quadratically  (or  at  least 
diminish)  as  the  flow  rate  is  increased. 

In  order  to  measure  the  ohmic  losses  across  the  cathode  and 
membrane,  simulations  were  performed  for  a  series  of  flow  rates 
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Fig.  2.  Comparisons  of  the  experimental  and  simulated  polarisation  curves  at  dif¬ 
ferent  electrolyte  inlet  velocities  and  different  POM  concentrations. 


between  50  and  1000  ml  min-1  with  a  POM  concentration  of  0.45 
M.  The  ohmic  losses  were  defined  as  IRohm/A  =  E  -\rj\-  VceU, 
where  bars  denote  spatial  averages  of  the  quantities  in  the  cathode 
and  A  is  the  active  surface  area.  For  a  sufficiently  high  flow  rate, 
mass  transport  losses  are  essentially  eliminated  and  Rohm  gives  a 
measure  of  the  overall  ohmic  resistance;  mass  transport  and  acti¬ 
vation  losses  in  the  conventional  anode  are  small  and  the  ohmic 
losses  are  incorporated  in  Ra.  The  change  in  Rohm  decreased  rapidly 
as  the  flow  rate  increased  from  50  and  250  ml  min-1 ,  beyond  which 
only  small  changes  were  observed.  A  flow  rate  of  1000  ml  min-1 
yielded  a  qualitative  prediction  of  Rohm  =  0.173  £2  cm2.  The  ohmic 
losses  in  a  conventional  cell  are  highly  dependent  on  the  mem¬ 
brane  conductivity  and  dimensions,  the  operating  temperature 
and  relative  humidity  [7],  the  volume  fraction  of  ionomer  in  the 
catalyst  layers  [14],  and  the  conductive  properties  and  dimen¬ 
sions  of  the  electrodes  and  plates.  Reported  values  for  the  ohmic 
resistance  in  state-of-the-art  conventional  cells  generally  fall  in 
the  range  50-1 00  m^  cm2  [13,14].  These  measurements,  taken 
from  current-interrupt  or  impedance  spectroscopy  experiments, 
cannot  be  compared  directly  to  the  Rohm  value  given  above, 
which  is  a  highly  qualitative  approximation.  Gasteiger  et  al.  [13] 
reported  a  cell  ohmic  resistance  of  45-55  m£2  cm2  (using  current- 
interrupt)  for  a  PEFC  employing  a  25  p,m  Nation  membrane  and 
operating  at  80  °C  with  fully  humidified  hydrogen  and  air  inlet 
streams.  In  their  cell,  which  used  a  microporous  layer  in  the  cath¬ 
ode,  mass  transport  limitations  were  not  dominant  even  at  the 
upper  limit  of  1 .2  A  cm-2 .  Makharia  et  al.  [  1 4]  estimated  the  ohmic 
resistance  of  a  PEFC  employing  a  50  p,m  thick  Nation  membrane  and 
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operating  at  80  °C  and  100%  relative  humidity  to  be  approximately 
60m£2cm2  (from  an  impedance  spectroscopy  analysis);  however, 
this  resistance  did  not  include  a  separately  calculated  catalyst-layer 
electrolyte  resistance  to  proton  transport,  which  was  estimated 
to  be  around  lOOm^cm2.  The  total  ‘ohmic’  resistance  of  around 
160  m£2  cm2  is  very  close  to  the  value  of  Rohm  given  above. 

The  condition  of  a  liquid-saturated  membrane  (at  the  cath¬ 
ode)  in  the  present  system  appears  to  maintain  a  high  membrane 
conductivity,  with  the  effect  of  Schroeder’s  paradox  likely  to  be 
a  contributing  factor.  Although  the  ohmic  losses  are  considerably 
reduced  in  a  conventional  cell  by  full  humidification  of  both  chan¬ 
nels,  severe  mass-transport  losses  caused  by  liquid-water  flooding 
can  limit  the  maximum  (limiting)  current  density  [7].  To  ameliorate 
this  problem,  micro-porous  layers  and  optimised  operating/flow 
conditions  are  used  [13,15].  Such  losses  are  eliminated  altogether 
in  the  present  cell  without  the  need  for  contact-angle  control  using 
hydrophobic  agents  and/or  additional  components  such  as  an  MPL 
(which  can  degrade  [16])  and/or  water  management  strategies,  as 
well  as  higher  flow  rates/pressures  in  the  channels.  As  is  seen  from 
the  simulated  curve  in  Fig.  2(b),  moderately  high  power  densities 
up  to  and  beyond  1.8  A  cm-2  are  readily  achievable  in  the  present 
unoptimised  cell. 

4.  Conclusions 

Experimental  and  modelling  results  for  a  novel  PEM  fuel  cell 
employing  a  solution  based  cathode  have  demonstrated  that  the 
cathode  overpotentials  under  typical  operating  conditions  can  be 
lower  than  those  in  a  cell  employing  a  standard  MEA.  Improve¬ 
ments  are  clearly  required,  nevertheless,  to  match  the  performance 
of  state-of-the-art  cells.  It  is  important  to  reiterate,  however,  that 
the  test  cell  studied  used  in  this  work  has  not  been  optimised  in 
terms  of  the  individual  components,  anode  humidification,  contact 
resistances,  POM  solution  composition  and  the  flow  distribution 
in  the  cathode,  amongst  other  considerations.  In  this  sense,  the 
results  are  highly  promising.  One  of  the  advantages  of  using  a  flow¬ 
through  electrode  is  the  ability  to  exert  control  over  the  reactant 
distributions  by  careful  design  of  the  flow.  In  conventional  cells, 
mass  transport  is  limited  by  the  diffusion  of  02  to  the  reaction  sites 
through  the  gas  diffusion  and  catalyst  layers,  which  is  more  difficult 
to  control  and  can  be  hindered  by  the  slow  gas-solid  reaction  and 


liquid-water  flooding  at  high  current  densities  and/or  high  relative 
humidities. 

The  model  results  in  this  paper  shows  that  the  losses  incurred  at 
the  anode  and  through  contact  resistances  (approximated  by  IRa ) 
are  non-negligible  and  that  the  effects  of  pEI  changes  on  the  reac¬ 
tion  kinetics,  Nernst  potential  and  DOR  are  potentially  significant. 
Work  is  ongoing  to  explicitly  incorporate  both  the  anode  and  the 
influence  of  pEI  and  to  investigate  the  stoichiometry  of  the  cath¬ 
ode  reaction  in  detail.  Future  studies  based  on  a  full-cell  model  will 
examine  voltage  losses  in  greater  detail  under  different  operating 
conditions  and  design  parameters. 
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